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INTRODUCTION 


There is great need for atmospheric trace constituent measurements with 
higher vertical resolution than attainable with passive radiometers. Infrared 
(IR) differential absorption lidar (DIAL) , which depends on Mie scattering from 
aerosols , has special advantages for tropospheric and lower stratospheric appli- 
cations and has great potential importance for measurements from Shuttle (ref. 1) 
and aircraft. DIAL data reduction involves comparing large amplitude signals 
which have small differences. The accuracy of the trace constituent concentra- 
tion inferred from DIAL measurements depends strongly on the errors in deter- 
mining the amplitude of the signals. Thus, the commonly used SNR expression 
(signal divided by noise in the absence of signal) is not adequate to describe 
DIAL-measurement accuracy and must be replaced by an expression which includes 
the random coherent (speckle) noise within the signal (refs. 2, 3, and 4). A 
comprehensive DIAL computer algorithm (ref. 5) is modified to include heterodyne 
detection and speckle noise. Results of a parametric study are presented and 
comparisons with direct detection are discussed. Examples are given for mon- 
itoring vertical distributions of O^, H^O, and NH using a ground-, aircraft-, 
or Shuttle-based pulsed tunable C0 2 laser DIAL system. 

ANALYSIS OF DIAL SENSITIVITY WITH HETERODYNE DETECTION 

The expectation value P of the number of measured photons from one coher- 
ence area of a scattering cell at range R of length Ar is 
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where r\ = detector efficiency, T - optical efficiency, E = laser energy, cT 
3 = 180° backscatter coefficient per length per steradian, Ar = cell length = — , 
T = integration time, V = frequency, A = transmitter area (= receiver area for 
heterodyne system with single detector) , 5 — extinction coefficient (total 
minus that of measured gas), G = absorption coefficient of measured gas, and 
p = density of gas being measured. 

In differential absorption, measurements are made of two frequencies 
selected to maximize signal and differential absorption of the species being 
measured while minimizing interference effects. The double ratio of signals at 
adjacent scattering cells at two frequencies yields information about absorp- 
tion in the region between the scattering cells. 
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where P^j = expectation value of measured return from cell j at frequency i, 
3j_j - backscatter coefficient from cell j at frequency i, Aa = absorption 
coefficient difference between frequencies 1 and 2 , A£ = extinction coeffi- 
cient difference between frequencies 1 and 2 (not including gas being measured) , 
i = 1 or 2 frequencies on or off absorption line, respectively, and j - 1 or 2 
for distances and R 2 , respectively. For scattering cells of equal length 
the resolution length R 2 - R^ equals the scattering cell length Ar. The term 
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In -5 q- — is a correction term due to changes in backscattering with fre- 
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quency across scattering cells. The term 2 / (A£) dr is a correction due to 

interferent species. Variation in backscatter and interferent absorption at 
the two frequencies will result in biases in the inferred concentrations. 

These biases can be reduced by careful frequency selection and partially cor- 
rected by using a priori information and auxiliary measurements. For an opti- 
cally thick species, a series of "on" frequencies is required to maximize sensi- 
tivity at various altitudes. 


The random error in p, which is calculated assuming that signal plus 
background are measured during each pulse and that background is measured 
between pulses and subtracted, will depend on the random error in the measure- 
ment of Pp-i- The uncertainty in the inferred concentration (5p ) 2 is given by 
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where N is the number of pulse pairs per measurement and SNR is the single 
pulse signal-to-noise ratio. For heterodyne detection, the major errors in 
P^j are due to quantum noise in the local oscillator and fluctuation noise in 
the return signal. 
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The heterodyne signal-to-noise ratio for a single coherence volume 

is limited to 1.0 due to the speckle noise in the return signal. The 
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number of coherence lengths per scattering cell is BT. The number of coher- 
ence areas viewed by the detector is M, where M is the ratio of receiver to 
transmitter area. For heterodyne detection an individual detector is required 
for each coherence area. The single pulse signal-to-noise ratio for determin- 
ing the signal from a scattering cell is then 
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where BTM is the number of statistically independent samples from a scat- 
tering cell for a single pulse. The post-detection bandwidth B is con- 
strained by matching with pulse duration (B < 1/Tp) and by the width of 
atmospheric spectral lines. The integration time T is constrained by verti- 
cal resolution requirements. Neglecting the system error terms, the random 
error in density for heterodyne detection with bandwidth B, integration time 
T, N pulse pairs, and M coherence areas with one detector per coherence 
area is 
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For N pulse pairs and M detectors, the SNR — — ■ -- - /BTNM is maximized 

for a given total laser energy per measurement (<* NMP-^j) when the laser energy 
per pulse per detector is selected so that P^j ~ BT. For that condition an 
approximate solution to equation (5) is 
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This equation illustrates that Sp is independent of p and proportional to 
1/Aa within the constraint that the two-way integrated absorption through the 
atmosphere does not limit Pj_j to less than BT. Equation (6) also indicates 
the strong effect of range resolution on measurement error. The number of 
pulse pairs required to maintain a constant measurement error is proportional 
to the inverse cube of the range resolution. 


COMMENTS ON DIRECT DETECTION 


Direct detection can have advantages over single-detector heterodyne 
detection when signal levels are high since direct detection allows averaging 
over multiple coherence areas with a single detector. The major disadvantage 
of direct detection is background noise, which is limited by optical filters 
(% lO 1 -^- Hz) compared with the electronic filters (10 6 to 10^ Hz) for heterodyne 
detection. Reducing the background noise by reducing the field of view will 
result in an increase in speckle noise. 


The direct detection signal-to-noise ratio can be written as a function 
of M, the number of coherence areas, and the total direct detected signal 
P^j = MPij. The SNR for direct detection (assuming zero detector noise, 
ref. 3) is 
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2 y. 

Equation (7) includes a speckle term p ij '/M, a Poisson term P^ jBT , and a 
background term (GDMT)BT. In equation (7) , G is the background signal in 
detected photons/sec/Hz per coherence area, M is the number of coherence 
areas, and D is the bandwidth in Hz of the optical filter. Maximization of 
the direct detection SNR with respect to the number of coherence areas occurs 
when M is selected so that the background and speckle terms are equal: 
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For this condition, the criterion for direct detection SNR to exceed the 
single-detector heterodyne large signal SNR of /bt is for P^. >2T/SBD. 


SENSITIVITY ANALYSIS FOR DIAL MEASUREMENT USING 
HETERODYNE DETECTION WITH A SINGLE DETECTOR 

A comprehensive computer algorithm is used to calculate the expectation 
values of . for various measurement conditions. Pressures, temperatures, 
and gas species densities are input from a midsummer midlatitude atmospheric 
model. Trace gas species densities can be modified using card inputs. Line 
absorption parameters are accessed from a comprehensive data base. Sources 
for line data are given in reference 6. At each altitude, molecular absorption 
at V]_ and V 2 is calculated for each species by summing contributions from 
absorption lines in the vicinity of the laser frequency. Lorentz , Voigt, or 
Doppler line shapes are used at appropriate altitudes. Water vapor continuum 
absorption is added to the line absorption. Extinction due to particulate and 
molecular scattering is summed with molecular absorption at each altitude to 
give the total loss in each scattering cell. The integrated two-way loss is 
calculated by summing contributions from altitude layers between the laser and 
the cell being considered. The backscattering coefficient 3 is calculated 
by combining a Rayleigh term, which is small at infrared frequencies, and a Mie 
term. Mie backscatter and extinction for the cases shown here are calculated 
using parameters of Deirmend j ian 1 s Haze L size distribution (ref. 7). The 
vertical aerosol distribution is based on reference 8 with a ground level 
concentration (350 particles/cm 3 ) corresponding to a 23-km visibility. 

Figure 1 shows the Mie contribution to the volume backscatter coefficient used 
in this study as a function of altitude for V = 927.61 cm - ' 1 '. 

For all simulations presented, the overall system efficiency, which is 
the product of optical and quantum efficiency, is set at 12.5%. The system 
errors due to amplification and digitization are set at 0.1%. For aircraft 
and ground-based cases, telescope area is 0.1 m^ and bandwidth 10 8 Hz; for 
Shuttle cases, telescope area is 1 m^ and bandwidth 10^ Hz. Bandwidths have 
been constrained by atmospheric line parameters and by laser limits. Pulse 
energies are selected within practical laser constraints to match signal to 
bandwidth and the number of pulse pairs N is adjusted to obtain reasonable 
measurement accuracy. The "on" and "off" frequencies have been selected to be 
within the range of a single rare isotope multiatmospheric CO laser line. Use 
of rare isotope lines minimizes interference by atmospheric C0 2 and use of 
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closely spaced "on" and "off" frequencies minimizes error due to variations in 
backscatter and interfering species absorption. 

Figure 2 displays O 3 concentration in STP-ppb and measurement error in 
the same units versus altitude for measurement from Shuttle at 250 km. The 
pulse energy is 5 joules per pulse with 1000 pulse pairs per measurement. 
Telescope area is 1 m^ , bandwidth is 10^ Hz, resolution is 1.5 km below 30 km 
and 3 km above 30 km. Simulations are for an "off" frequency of 1058.01 cm ” 1 
and "on" frequencies of 1058.17, 1058.11, and 1058.19 cm -1 . It should be 
noted that a fine tuning of the "on" frequency results in a sharp variation 
of altitude at which the best measurement can be accomplished. 

Figure 3 displays ozone concentration and error in STP-ppb for a measure- 
ment from aircraft at 10.5 km looking either upward or downward. Pulse energy 
is 0.05 joules, telescope area is 0.1 m 2 , bandwidth is 10 s Hz, and "on" fre- 
quency is 1058.20 cm . in the case of aircraft measurements the "on" frequency 
selection is simplified since the signal does not have to pass through the ozone 
bulge. Two cases are presented in order to illustrate the strong influence of 
range resolution on number of pulses required for a measurement. The solid er- 
ror line is for 3000 pulse pairs with a resolution of 0.5 km below 6 km, 1 km 
from 6 to 10 km, and 3 km above 10 km. The dashed error line is for 100 pulse 
pairs with a resolution of 1.5 km. The errors are nearly identical at low al- 
titude, which is expected (see equation ( 6 )) since the number of pulses has been 
increased by approximately the cube of the inverse of the range resolution. 

Figure 4 displays NH 3 concentration and error in STP-ppb for 300 pulse 
pairs with operation from Shuttle at an altitude of 250 km with 1 joule per 
pulse and from aircraft at an altitude of 10 km with 0.05 joules per pulse. 

Since NH 3 is not heavily attenuated, only one set of frequencies is required. 

The resolution for both cases is 1.5 km. 


Figure 5 displays NH 3 concentration and error in STP-ppb for a ground- 
based measurement. The vertical NH 3 distribution is the same as shown in 
figure 4. For this case, the telescope area is 0.1 m 2 , bandwidth is 10 8 Hz, 
"on" frequency is 927.32 cm -1 , and "off" frequency is 927.61 cm -1 . The error 
curves are for a range resolution of 1.5 km using 300 pulse pairs. The two 
error curves are for laser energies of 0.05 J and 0.5_J. Note that for the 

cases in which E = 0.5 J, the large signal limit 


p ij + 


BT 


is achieved 


throughout the entire measurement range and therefore the error is nearly 
independent of altitude. 


Figure 6 is for the same conditions as shown in figure 5 with the excep- 
tion that the NH 3 concentration has been reduced by a factor of 5. The error 
is plotted only for the E = 0.05 J case. Comparison of figure 6 with the 
0.05 joule curve in figure 5 illustrates the effect of concentration on 
measurement error. At low altitude, P^j >BT and the measurement errors are 
independent of concentration. At higher altitudes P-^j <BT and the increased 
integrated absorption at the higher concentration (fig. 5) results in an 
increased measurement error for that case. These NH 3 distributions are 


5 



typical values from Langley Research Center ground-based infrared heterodyne 
radiometer measurements (ref. 9) . 


Figure 7 displays percent error in measurement of water vapor concentra- 
tion versus altitude for measurement from 250 km for a midsummer, midlatitude 
water vapor distribution. Measurement conditions are 1 joule per pulse, 100 
pulse pairs, 10* 7 Hz bandwidth, and 1 m 2 telescope. Two "on" frequencies are 
used at 948.25 and 948.30 cm” 1 . The "off" frequency is at 948.35 cm” 1 . Trop- 
ospheric water vapor is particularly amenable to measurement from space since 
the concentration rapidly increases towards the ground. Thus, absorption is 
high at low altitudes and the integrated absorption from space to low alti- 
tudes is relatively low. For measurement of water vapor at higher altitudes, 
the number of pulses would have to be increased. 


CONCLUDING REMARKS 

High vertical resolution measurement of atmospheric trace species can be 
achieved using CO^ laser DIAL with heterodyne detection. This study indicates 
that maximum sensitivity at minimum laser energy per measurement requires 
multiple pulse operation with the energy per pulse selected so that the measured 
photon rate is approximately equal to the detector IF bandwidth. Measurement 
sensitivities can be maximized and interference effects minimized by fine ad- 
justment of measurement frequencies using the tunability of high-pressure lasers. 
The use of rare isotope lasers minimizes loss due to CO^ atmospheric absorption. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 5, 1980 


6 



REFERENCES 


■i 

1 


1. Shuttle Atmospheric Lidar Research Program - Final Report of Atmospheric 

Lidar Working Group. NASA SP-433, 1979. 

2. Jakeman, E. ; Oliver, C. J. ? and Pike, E. R. : Optical Homodyne Detection. 

Adv. Phys., vol. 24, no. 3, May 1975, pp. 349-405. 

3. Elbaum, Marek? and Teich, Malvin Carl: Heterodyne Detection of Random 

Gaussian Signals in the Optical and Infrared: Optimization of Pulse 

Duration. Opt. Commun. , vol. 27, no. 2, Nov. 1978, pp. 257-261. 

4. Rye, B. J. : Differential Absorption Lidar System Sensitivity With 

Heterodyne Reception. Appl. Opt., vol. 17, no. 24, Dec. 15, 1978, 
pp. 3862-3864. 

5. Remsberg, Ellis E. ; and Gordley, Larry L. : Analysis of Differential 

Absorption Lidar From the Space Shuttle. Appl. Opt., vol. 17, no. 4, 

Feb. 1978, pp. 624-630. 

6. Park, J. H. : Optical Measurement in the Middle Atmosphere. Pure & Appl. 

Geophys., vol. 117, no. 3, 1978-1979, pp. 395-429. 

7. Deirmend jian, D. : Electromagnetic Scattering on Spherical Polydispersions. 

American Elsevier Publishing Co., New York, NY, 1969. 

8. McClatchey , Robert A.; and Selby, John E. A.: Atmospheric Attentuation of 

Laser Radiation From 0.76 to 31.25 yra. AFCRL-TR-74-0003 , U.S. Air Force, 
Jan. 1974. 

9. Hoell, James M. ? Harward, Charles N. ; and Williams, Burnie S. : Remote 

Infrared Heterodyne Radiometer Measurements of Atmospheric Ammonia Pro- 
files. Geophys. Res. Lett., vol. 7, no. 5, May 1980, pp. 313-316. 




7 




Figure 1.- Aerosol volume backscattering coefficient as a function of 

altitude for V = 927.61 cm“-*-. 



Ar = 1.5 km Ar = 3 km 


Figure 2.- Ozone measurement from 250 km. 
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ALTITUDE IN KILOMETERS 


Figure 3.- Ozone measurement from 10.5 km. 



Ar = 1. 5 km 

Figure 4.- NH^ measurement from 10 and 250 km. 
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ALTITUDE IN KILOMETERS 
Ar = 1.5 km 

Figure 5.- NH^ measurement from ground for E = 0.5 and 0.05 J 



ALTITUDE IN KILOMETERS 


Figure 6.- NH 0 measurement from ground for E = 0.05 J. 
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Ar = 1. 5 km 


Figure 7.- Water vapor measurement from 250 1cm. 
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